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The insulator to metal transition in LaTiOs induced by La substitution via Sr is studied within 
multi-band exact diagonalization dynamical mean field theory at finite temperatures. It is shown 
that weak hole doping triggers a large interorbital charge transfer, with simultaneous electron and 
hole doping of t2g subbands. The transition is first-order and exhibits phase separation between 
insulator and metal. In the metallic phase, subband compressibilities become very large and have 
opposite signs. Electron doping gives rise to an interorbital charge flow in the same direction as 
hole doping. These results can be understood in terms of a strong orbital depolarization. 
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The remarkable sensitivity of the electronic and mag- 
netic properties of transition metal oxides to small 
changes of parameters such as temperature, pressure, 
or impurity concentration has attracted wide attention 
during recent years A phenomenon of particular 

interest is the filling-controlled metal insulator transi- 
tion in the vicinity of integer occupancies of partially 
filled valence bands. For instance, LaTiOa is known 
to be a Mott insulator which becomes metallic if La 
is replaced by a few percent of Sr ions [3, 0]. This 
kind of doping-driven insulator to metal transition has 
been discussed theoretically in terms of single-band mod- 
els 0, i, S 0, i, i, M, El, [IljJIlor multi-band models 
based on degenerate orbitals \l4 . TEl IT6| . To our knowl- 
edge, only one calculation for non-identical subbands ex- 
ists so far p7j . 

It has recently been shown that orbital degrees of free- 
dom play a crucial role in the Mott transition of materials 
consisting of non-equivalent partially occupied subbands. 
For example, whereas cubic SrVOa {3d^ configuration) is 
metallic, iso-electronic LaTiOa exhibits a Mott transition 
because non-cubic structural distortions lift the degener- 
acy among t2g subbands ).18j . Local Coulomb interac- 
tions strongly enhance this orbital polarization, so that 
the electronic structure of the insulating phase consists 
of a single nearly half-filled band split into lower and 
upper Hubbard peaks, and two nearly empty t2g sub- 
bands. Similarly, the monoclinic symmetry of the Mott 
insulator V2O3 (Sd^) yields a pair of nearly half-filled 
subbands (with lower and u Tppe r Hubbard peaks) and a 
single nearly empty t2g band [l9l.[20|. Also, the insulating 
phase of the layered perovskite Ca2Ru04 {M'^) consists 
of a fully occupied t2g band of predominant dxy character 
and two half-filled dj^z^y^-like bands [2l|. In view of the 
nearly complete orbital polarization of the Mott phase 
of these systems it is evidently necessary to go beyond 
single-band or degenerate multi-band models to describe 
the effect of electron or hole doping. 

In the present paper we investigate the destruction of 
the Mott phase of LaTiOs due to replacement of small 
amounts of La by Sr. To account for the orbital de- 
grees of freedom during hole doping we make use of a 



recent multi-band extension of dynamical mean field the- 
ory based on exact diagonalization (2^ . This approach 
has been shown to be highly useful for the study of strong 
correlations in several transition metal oxides [21I [2^. ^2^^] . 
Since it does not suffer from sign problems the full Hund 
exchange including spin-flip and pair-exchange contribu- 
tions can be taken into account. Also, relatively large 
Coulomb energies and low temperatures can be reached. 

Naively, one might expect that all subbands are af- 
fected similarly by the doping process, perhaps weighted 
by their relative occupancies. The result of this work is a 
completely different picture. As we show below, the addi- 
tion of only a few percent of holes in the t2g bands triggers 
a substantially larger interorbital charge rearrangement, 
with a strong flow of electrons from the nearly half-filled 
subband towards the nearly empty ones. In other words, 
weak "external" hole doping of the t2g bands causes a 
much larger "internal" simultaneous electron and hole 
doping of different subbands. Thus, the density-driven 
insulator to metal transition in Lai-ajSr^jTiOa is accom- 
panied by a significant reduction of orbital polarization. 
We also show that electron doping of the t2g bands, which 
could be achieved by substitution of La by Zr, leads to 
a strong charge flow in the same direction as hole dop- 
ing, i.e., from the nearly half-filled t2g band to the nearly 
empty ones. 

As shown by Pavarini et al. [18], the octahedral distor- 
tions due to the orthorhombic lattice of LaTiOs ensure 
that the t2g density of states is not diagonal within the 
xy,xz,yz basis. It is nonetheless possible to construct 
from the <2g orbitals a new basis which is nearly diago- 
nal [l^l. We denote this basis here as Ug and e'g. The 
corresponding densities of states obtained from ab ini- 
tio calculations within the local density approximation 
(LDA) are shown in Fig. 1. The centroid of the Ug den- 
sity lies about 200 mcV below that of the two nearly 
degenerate e' contributions. Accordingly, the subband 
occupancies (per spin band) in the absence of correlation 
effects are Ua ~ 0.23 and « 0.135. This t2g crys- 
tal field splitting compares well with recent experimental 
estimates based on X-ray studies [2^ [2^ . 

To account for local Coulomb interactions among the 
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FIG. 1: (Color online) Density of states components of 
LaTiOa in nearly diagonal t2g subband representation [2^ . 
Solid (red) curve: Ug states, dashed (blue and green) curves; 
e'„ states; Ef ~ 0. 
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FIG. 2: (Color online) Subband occupancies (per spin band) 
of Lai-ajSr^iTiOa as functions of Coulomb energy. Solid (red) 
curves: 3d^ {x — 0); dashed (blue) curves: 3d°'^^ {x = 0.05); 
long-dashed (green) curves: 3dP'^° (x = 0.1). 
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Ti "id electrons we use finite temperature dynamical mean 
field theory (DMFT) [27, 28, 29, 30] combined with exact 
diagonalization (ED) [3l|. It was recently shown [22] that 
this method can be generalized to multi-band materials 
by computing only those excited states of the impurity 
Hamiltonian that are within a narrow range above the 
ground state, where the Boltzmann factor provides a suit- 
able convergence criterion. In the present case, each of 
the three i2g impurity orbitals couples to two bath levels, 
giving a cluster size rig = 9 with maximum Hamiltonian 
dimension 15.876. Exploiting the sparseness of the im- 
purity Hamiltonian, a limited number of excited states 
can be computed very efficiently by using the Arnoldi 
algorithm [32| . Since local Coulomb and exchange inter- 
actions couple the baths of the three t^g orbitals, the level 
spacing between excited states is very small so that finite 
size effects are much smaller than for a single impurity 
level coupled to only two bath levels. 

Fig. 2 shows the variation of the ag and sub- 
band occupancies with Coulomb energy. Full Hund ex- 
change is included, and the exchange energy is held fixed 
at J = 0.65 eV [l^. The inter-orbital Coulomb en- 
ergy is given by J7' = \J — 2J . The temperature is 
T = 10 ... 20 meV. For simplicity, the small difference 
between the e'g densities is neglected and their average 
density is used instead. Integer and non-integer occu- 
pancy are seen to lead to qualitatively different behavior. 
For n ~ 2na + 4ne = 1, local Coulomb interactions grad- 
ually suppress orbital fluctuations up to f/ w 4.5 eV. In 
the subsequent narrow range up to ?7 = 5.0 eV, the Ug 
band is rapidly becoming nearly half-filled and the two 
e'g bands nearly empty. This result supports the picture 
presented in Refs. [l8|, [IJ, where for U = 5.0 eV, Ising 
exchange, and T = 100 meV, LaTiOa was found to be a 
Mott insulator, with nearly empty e'g subbands, and the 
half-filled ag band split into lower and upper Hubbard 
peaks [33} . 

For hole doping with n = 0.95 and 0.90, on the other 



hand, the orbital polarization increases smoothly with in- 
creasing Coulomb energy, without any sign of the critical 
behavior that is characteristic of the pure La compound. 
Evidently, the rapid suppression of orbital fluctuations 
near U — 5.0 eV is a phenomenon associated with the 
integer occupancy of the t2g bands. To illustrate the de- 
struction of the Mott phase of LaTiOs via doping with 
Sr we assume here that the primary effect is caused by 
the change of band filling. Additional modifications of 
the Ug and e'g density of states due to the slightly smaller 
size of Sr compared to La are neglected. 

The unexpected and striking feature of these results 
is that the charge rearrangement among t2g orbitals is 
much larger than the number of holes induced via La — ^ 
Sr substitution. For instance, for x = 0.05 (n = 0.95) the 
total ag occupancy (both spins) changes from near unity 
to 0.64, whereas the total e'g occupancy (four spin bands) 
changes from near zero to 0.31. Thus, the internal charge 
transfer from e'g to ag bands is six times larger than the 
external charge transfer due to doping with Sr. Clearly, 
the destruction of the Mott phase of LaTiOs via hole 
doping does not proceed via approximately equal partic- 
ipation of all subbands. Instead, weak external hole dop- 
ing generates a much larger simultaneous electron and 
hole doping in different t2g subbands. 

Since the density of states of LaTiOa is not symmetric 
and the t2g occupation is far from one-half, it is inter- 
esting to consider also the case of electron doping, for 
example, via substituting La by Zr. Because the ag band 
is nearly full, one might expect the donated electrons to 
populate the nearly empty e'g bands, with little other 
changes. As in the hole doping case, we find an entirely 
different picture. For example, for 0.10 electron doping, 
0.26 electrons are injected into the e'g bands and 0.16 
holes into the ag band. Thus, the actual number of elec- 
trons transfered to the e'g bands is 2.6 times larger than 
externally supplied. Of course, this can only be achieved 
via simultaneous hole doping of the nearly filled ag band. 
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FIG. 3: (Color online) Red and blue dots: doping induced 
change of total ag and e'g subband occupancies, Aa and Ae, 
respectively, for hole doping (c ~ —0.05) and electron doping 
(c ~ 0.1); U — 5.0 eV. In both cases, holes are injected into 
the nearly half-filled ag band (red arrows) and electrons into 
the nearly empty e'g bands (blue arrows). The green arrows 
denote the net doping. 
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FIG. 4: (Color online) Solid (red) and dashed (blue) curves: 
ag and e!g subband occupancies (per spin band) of LaTiOa 
as functions of chemical potential for U — 5.0 eV. Long- 
dashed (green) curve: average occupancy Uav = (na -f 2ne)/3. 
fi < 1.5 eV corresponds to hole doping, /i > 2.5 eV to electron 
doping. The arrows mark the hysteresis behavior for increas- 
ing and decreasing n, indicating the first-order nature of the 
density-driven Mott transition. 



These results are illustrated in Fig. 3 which shows the 
total electronic charge injected into the e'g bands (Ag — 
Aue > 0) and the total hole charge injected into the ag 
bands (A^ = 2na — 1 < 0) for net doping concentrations 
c w —0.05 and c w 0.1, where c = A^, + = n—1. Both 
electron and hole doping are seen to cause a strong orbital 
depolarization, implying an internal charge flow between 
t2g bands that is significantly larger than the externally 
supplied number of electrons or holes. Note also the lack 
of symmetry between electron and hole doping. 

Fig. 4 shows in more detail the charge transfer among 
t2g bands as a function of chemical potential. The 
Coulomb and exchange energies are held constant at 
values appropriate for LaTiOa. Hole doping occurs for 
/i < 1.5 eV, electron doping for > 2.5 eV. Because of 
the insulating gap, the charge compressibility k — dn/dfi 
vanishes in the intermediate range, i.e., the chemical po- 
tential can be varied without changing the subband occu- 
pancies. The small deviations from Ua = 0.5 and Ug = 0.0 
in this region presumably are due to the finite tempera- 
ture and ED finite size effects. The hysteresis loops near 

— 1.5 eV and ^ — 2.5 eV demonstrate that the hole 
and electron driven Mott transitions at finite T are first 
order. Phase separation, where an insulating solution 
with n — 1 coexists with a metallic solution for n ^ 1, 
is seen to occur at both transitions. Moreover, the com- 
pressibility gets very large as the upper and lower critical 
values of fj, at both loops are approached. Note that the 
hysteresis behavior of the ag and e'g subbands is much 
more pronounced than that of the average t2g charge. 
Accordingly, the subband compressibilities Ki — drii/dfj, 
have opposite signs and their absolute values are much 
larger than the average t2g compressibility. These re- 
sults show that the orbital degrees of freedom lead to a 
non-trivial generalization of the one-band picture close to 
half-filling g, H, i, 0, i, i, 0, [m, E [il and of the multi- 



band picture based on equivalent orbitals [TJ, [T^ [lB| ■ 

Finally, to make contact with the photoemission data 
for Lai-^Sr^TiOa [H, we show in Fig. 5 the quasi-particle 
spectra for ti = 1 (x = 0) and n = 0.95 {x = 0.05). 
For simplicity, we give here the ED cluster spectra which 
can be evaluated directly at real w, without requiring ex- 
trapolation from Matsubara frequencies. At integer oc- 
cupancy the system is insulating. The excitation gap is 
formed between the lower Hubbard peak of the a„ band 
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FIG. 5: (Color online) Quasi-particle spectra of LaTiOa (up- 
per panel) and Lao.gsSro.osTiOs (lower panel) for U — 5 eV. 
Solid (red) curves: Ug states, dashed (blue) curves: e'g states; 
Ef = 0. 
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and the empty bands, in agreement with Ref. [23 |. 
The hole doped system, on the other hand, is metalhc, 
with conduction states stemming from both ag and 
bands. If we associate the peak near Ep with coherent 
states, the large peak at 1.1 eV binding energy, i.e., below 
the bottom of the LDA density of states, represents the 
incoherent weight related to the lower Hubbard band. It 
has mainly ag character, with weak e'g admixture. The 
position of this feature is in good agreement with the 
photoemission spectra for Lao.gsSro.osTiOs 0]. More- 
over, the ratio between incoherent and coherent weight 
is roughly 3:1. This is also consistent with the photoe- 
mission data which reveal substantially more incoherent 
than coherent emission. 

In summary, we have shown that multi-band ED 
DMFT is a powerful method to study the subtle in- 
terorbital charge rearrangement triggered by the density- 
driven Mott transition in Lai-^rSr^jTiOs. In particular. 



we have shown that doping with Sr gives rise to a first- 
order insulator to metal transition that is accompanied 
by a significant orbital depolarization. Injection of only a 
few percent of holes into the Ti t2g bands leads to a much 
larger simultaneous electron and hole doping of different 
subbands. The transition exhibits separation between 
insulating and metallic phases. In the latter, subband 
charge compressibilities are very large and have opposite 
signs. Moreover, electron doping generates an internal 
charge flow in the same direction as hole doping. Since 
other Mott insulators such as V2O3 and Ca2_a;Sra;Ru04 
also exhibit nearly complete orbital polarization, it is 
likely that doping these systems with electron or hole 
donors will lead to a similar enhancement of orbital fluc- 
tuations as discussed here for Lai_a;Sr^Ti03. 

I like to thank Eva Pavarini for the density of states 
distributions shown in Fig. 1, and Theo Costi for useful 
discussions. 
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